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Abstract. A microscopic theory for the interaction of carriers with LO phonons is used to study the 
ultrafast carrier dynamics in nitride-based semiconductor quantum dots. It is shown that the efficiency 
of scattering processes is directly linked to quasi-particle renormalizations. The electronic states of the 
interacting system are strongly modified by the combined influence of quantum confinement and polar 
coupling. Inherent electrostatic fields, typical for InGaN/GaN quantum dots, do not limit the fast scattering 
channels. 
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1 Introduction 

In the rapidly evolving field of semiconductor quantum 
dots (QDs), which provide many new applications in op- 
toelectronics, nitride-based systems are of growing interest 
due to their large range of possible emission frequencies.^ 
I2B31 The interaction of carriers with LO phonons is an 
important scattering process which determines the redis- 
tribution of carriers between the localized QD states and 
between QD and wetting-layer (WL) states. Furthermore, 
these interaction processes lead to optical dephasing which 
is of central importance in coherent optical experiments 
and determines the homogeneous broadening of absorp- 
tion and gain spectra. In nitride-based QDs, the role of 
interaction with LO phonons is increased due to the inter- 
mediate polar coupling strength in this material system. 
On the other hand, for QDs with wurtzite crystal structure 
grown along the c-axis, strong built-in fields modify the 
electronic states via the quantum-confined Stark effect. 0| 
The aim of this paper is to determine the efficiency of 
carrier-phonon interaction under these conditions. 

While a simplified treatment of carrier-phonon inter- 
action in QDs, based on Fermi's golden rule, leads to the 
prediction of a phonon-bottleneck (see |5] and references 
therein), recently it has been emphasized that perturba- 
tion theory is not applicable for QD systems and that the 
interaction needs to be treated in the polaron picture. 
[7] A quantum-kinetic theory has been used to describe 
carrier scattering processes due to QD polarons.JB] We 
extend this model by including also the optical genera- 
tion of carriers by means of a short laser pulse. Naturally 
linked to this is the dephasing of the coherently driven 
interband polarization, which is described on the same 
quantum-kinetic level as the carrier scattering. Further- 



more, we study intermediate polar coupling and include a 
self-consistent description of the quantum-confined Stark 
effect for nitride-based QDs. 



2 Quantum-dot model 

We consider a system of self-assembled InGaN/GaN QDs, 
grown on the c-plane in wurtzite crystal structure where 
the spontaneous and piezoelectric polarization along the 
c-axis cause large electrostatic fields. The calculation of 
single-particle properties for these QDs is a challenging 
task on its own. Considerable recent progress has been 
made in k ■ p and tight-binding models for realistic QD 
geometries of various material systems to determine the 
single-particle wave functions and energies and the result- 
ing free-particle optical transitions. |51ll0llllllT^] With the 
investigation of many-body interaction effects for given 
single-particle states, this paper follows a different aim. 
Previous calculations of carrier-carrier [T5ll 1 4j and carrier- 
phonon[Hj scattering in QDs showed that particular details 
of the wavefunctions are less important, while results are 
more sensitive to the energy stucture (level spacing and 
degeneracy) . 

For our purpose, we use a simplified model for the 
single-particle states of the QD system based on a sep- 
aration ansatz 3>(r) = f(p) £(z) for the wave functions. 
The in-plane part <p{p) describes either the localized QD 
states or the delocalized WL states which form a con- 
tinuum at higher energies. In growth direction, perpen- 
dicular to the WL plane, the carrier confinement poten- 
tial is modified by intrinsic electrostatic fields. Due to the 
quantum-confined Stark effect, charges become partly sep- 
arated and a screening field builds up which reduces the 
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Fig. 1. Confinement potential U(z) (solid line) and wave func- 
tion |£(z)| 2 (dashed line) for electrons (a) and holes (b) in 
growth direction for an electrostatic field of 2 MV/cm. 



intrinsic electrostatic field. Therefore £(z) has to be cal- 
culated from a self-consistent solution of the Schrodinger 
and Poisson equation. For details on the QD model, the 
calculation of wavefunctions and material parameters, see 
d 

We consider QDs with two confined shells which due to 
their angular momentum properties are referred to as the 
s-shell (ground state) and the p-shell (excited state). The 
level spacing between s- and p-states is 1.2 hui^o for elec- 
trons and 0.4 /zwlo for holes. For the intrinsic electrostatic 
field, we consider a typical field strength of 2 MV/cm. 



3 Many-body theory 

A perturbative treatment of the carrier interaction with 
LO phonons is only possible when modifications of the 
single-particle states and energies due to the polar cou- 
pling remain small. In QDs the localized states with dis- 
crete energies are subject to hybridization effects even in 
materials with weak polar coupling. 6,7,18] The QD po- 
laron states are broadened due to the interaction with car- 
riers in the WL continuum, which ensures dissipation. [S] 
Quasi-particle renormalizations can be described using the 
retarded Green's function, which obeys the Dyson equa- 
tion 



C% r (t) = S(r) 



+ f dr' £ G c /(r - r>) D%> (r - r') Q R (r') , U) 

JO a 







for electrons in the conduction band and a similar equation 
with _D VV,< for the valence band. This equation describes 
the polaron in the self-consistent, diagonal, random-phase 
approximation (RPA) |15ll5] . The single-particle energies 
for either localized or extended states are given by , 
with A = c,v being the band index. Here, the correspond- 
ing wave functions (which incorporate effects of the built- 
in field) enter via the interaction matrix elements *. 



They are included in the phonon propagator 
i 

x [ A^ LO e ±iWLOT + (1 + JV LO ) e^° T ] , (2) 

where A^o is a Bose-Einstein function for the population 
of the phonon modes (assumed to act as a bath in thermal 
equilibrium) and cjlo is the LO-phonon frequency. The 
Fourier transform of the retarded GF directly provides the 
spectral function G a (w) = —2 Im G„(w) which reflects the 
density of states (DOS) for a carrier in state a. 

For non-interacting carriers, the spectral function con- 
tains 5-functions at the free-particle energies, reflecting 
the infinite lifetime of these states (if recombination on 
a much longer timescale is neglected). The spectral func- 
tions for carriers in the QD-WL system interacting with 
lattice distortions are shown in Fig. For a WL electron 
with zero momentum (at the band edge), the 5-like spec- 
tral function in the absence of interaction (indicated by 
a vertical line in Fig. Ek) is broadened. The correspond- 
ing finite lifetime of the state is due to possible scattering- 
events which are self-consistently calculated in terms of 
renormalized states and in arbitrary order in the interac- 
tion by means of a direct numerical solution of Eqs. (1) 
and (2). Furthermore, the energy lowering of the interact- 
ing state (polaron shift) and weak LO phonon satellites 
are obtained. 

For the discrete QD states, Figs. Et") and c) show that 
the free-particle picture is strongly modified. On this level, 
a direct comparison between the the weak-polar coupling 
material system InGaAs (Fig. 1 of Ref. [5]) and the In- 
GaN system with intermediate polar coupling is possi- 
ble. In the first case, one can clearly distinguish between 
phonon replicas of the QD states which are spaced by 
the LO-phonon energy hu^o around the central peak, and 
the hybridization effect which leads to a splitting of these 
peaks that is smaller than the LO-phonon energy (when 
the detuning between the QD level spacing and the LO- 
phonon energy is small). In the InGaN material system, 
the hybridization is much stronger and one can no longer 
distinguish between phonon replicas and their splitting. 
Moreover, the obtained rich multi-peak structure provides 
a stronger overlap between the spectral functions which 
increases the scattering efficiency discussed below. Broad- 
ening of the resonances is provided by the coupling to the 
WL states and the finite LO-phonon lifetime (assumed to 
be 5 ps). 

In the markovian limit, the transition rate between 
two states a and /3 is proportional to the overlap A a p = 
J dio G a (u) G p{oj ±wlo) • In the free-particle picture, con- 
sidering the 5-like DOS, the results of Fermi's golden rule 
are recovered. The level-spacing has to match exactly the 
LO-phonon energy Swloj otherwise the overlap vanishes 
and carrier transitions are not possible. In the polaron 
picture, strong overlap between the spectral functions of 
the QD states as well as between the QD and WL states 
is obtained which already indicates the possibility of ef- 
ficient scattering processes. On ultrafast timescales, the 
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Fig. 2. Spectral function for a WL electron with zero momen- 
tum (a) and for a QD electron in the p-shell (b) and in the s- 
shell (c) . Vertical lines indicate the positions of the 5-functions 
for a non-interacting carriers. 



Markov approximation is not valid and quantum-kinetic 
effects additionally increase the scattering rates. 

To study the carrier relaxation from p-shell into s-shell, 
we consider the optical excitation of carriers with a short 
(100 fs) laser pulse having spectral overlap only with the p- 
shell. The coupled carrier and polarization dynamics are 
described by the kinetic equation for the single-particle 
density matrix p a (t), which contains optically induced in- 
terband transition amplitudes ip a and occupation proba- 
bilities /£' v according to ^B] 



d_ 

dt 



d_ 

di 



11 

r« fi 



fir 



-m 



dt 



Pa(t) 



coll. 



(3) 



In this matrix notation of the semiconductor Bloch equa- 
tions, Hartree-Fock renormalizations due to the carrier 
Coulomb interaction arc included in the single-particle en- 
ergies = e x + J2 
Q a = dE{t) + £ 



XX (fa ~ ^Av) and in the Rabi energy 
■ a pip a - Here d is the dipole coupling, 



f3 V al 3 

v 



E(t) the laser pulse and are the matrix elements of 
the Coulomb interaction :1~4*|. The quantum kinetic de- 
scription of the carrier-phonon interaction, based on the 



generalized Kadanoff-Baym ansatz JJj] , 



dt 



- [Sx,-p X ^(t')}pf(t')D^>(t,t') 
+ pfW) Kx>-pf(t')\ D^<(t,t>) 
+ G x a R (t,t') [<' R (M') 

+ [Sx,-p X nt')}p$ X '(t')D^^(t,t') 
- (frtf) [S„y-P$ x '(t')} f)]} - 

(4) 

represents a generalization of the formulation in Ref. [S] 
to the coupled carrier and polarization dynamics and de- 
scribes scattering, dephasing, and energy renormalization 
contributions. The RPA self-energy was used, and Hartree 
terms were neglected, since we consider only low-density 
populations. In the markovian limit, the slowly varying 
part of p(t') is taken out of the integral and memory ker- 
nels of the type A% = J dt'G x a R {t') [G^ K {t')} * e ± ^ t ' 
determine the scattering efficiency between two states a 
and (3 according to the overlap of spectral functions in the 
frequency domain, as stated above. In the non-markovian 
case, the memory kernel contains also the single particle 
density matrix. This gives rise to memory effects since the 
actual time evolution depends explicitly on the history of 
the system. 

The considered level-spacing of 1.2 ftw^o for electrons 
and 0.4 fi^LO f° r boles provides no carrier transitions in a 
free-particle picture. Results of the quantum-kinetic calcu- 
lation are shown in Fig-El Due to pulse excitation, first the 
p-shell is populated. The carrier relaxation is so fast that 
already during the optical pulse a strong redistribution 
from the p-shell to the s-shell takes place. The ultrafast 
scattering processes are connected with strong dephasing 
of the coherent interband polarization, which is included 
on the same quantum-kinetic level in Eq. (4). The strong 
dephasing leads to a time evolution of the interband po- 
larization, which practically follows the optical pulse (not 
shown). Correspondingly, the coherent regime is limited to 
the duration of the 100 fs pump pulse. A transient strong 
increase of the total carrier density is obtained during the 
coherent regime, as can be seen in Fig. 3 b). 

The time-domain oscillations in Fig. 3a) indicate the 
strong coupling between carriers and phonons. The oscil- 
lation frequency is connected to the hybridization in the 
DOS. The phenomenon is analogous to Rabi-oscillations, 
though not due to coupling of carriers to the optical field, 
as recently discussed for QDs in Refs. |17U18| . but rather 
due to coupling of carriers to the bosonic field of the 
phonons. Note that the oscillations appear after the op- 
tical pump pulse in the incoherent regime. Hence they 
are not connected to optical coherence. At larger times, 
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Fig. 3. (a) Time evolution of the electron population in the 
p-shell (dashed line) and s-shell (solid line) after optical exci- 
tation with a 100 fs pulse resonant with the p-shell. (b) Total 
carrier density (solid line), and carrier density in the WL (dot- 
ted line) and QD (dashed line) versus time. 

the possibility of carrier transitions between QD and WL 
states leads to a redistribution of carriers from the QD 
states into the WL. The corresponding reduction of the 
QD population in connection with the increasing WL car- 
rier density can be seen in Figs. 3 a) and b), respectively. 
Of course, the total carrier density remains constant after 
the pulse. 

The ultrafast scattering and dephasing processes are 
the direct result of the polaronic renormalizations. Only 
within a non-markovian calculation the above mentioned 
oscillations in connection with the strong-coupling regime 
for the carrier-phonon interaction can be obtained. [5] In 
this respect, there is no qualitative difference between the 
InGaAs and InGaN material systems. In the latter case, 
the larger polar coupling increases the scattering efficiency 
while the electron and hole charge separation induced by 
the built-in fields reduces only inter-band but not intra- 
band processes. 

In the presented example, a weak excitation pulse has 
been used to justify the restriction of the interaction to LO 
phonons. At elevated excitation densities, carrier-carrier 
scattering additionally contributes to fast scattering processes. HH 

4 Conclusion 



carriers with LO phonons. Therefore, Fermi's golden rule 
is not applicable and a quantum-kinetic description is nec- 
essary, which provides ultrafast scattering channels. 
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It has been shown that in nitride-based QD systems with 
intermediate polar coupling strong modifications of the 
free-particle states are obtained due to the interaction of 



